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We investigated by numerical simulations the effects of Schmidt number on passive scalar transport 
in forced compressible turbulence. The range of Schmidt number (Sc) was 1/25 ~ 25. In the 
inertial-convective range the scalar spectrum was seemed to obey the k ~ 5 ' 3 power law. For Sc S> 1, 
there appeared a fc” 1 power law in the viscous-convective range, while for Sc -C 1, a fc” 17 / 3 power 
law was identified in the inertial-diffusive range. The scaling constant computed by the mixed 
third-order structure function of velocity-scalar increment showed that it grew over Sc, and the 
effect of compressibility made it smaller than the 4/3 value from incompressible turbulence. At 
small amplitudes, the probability distribution function (PDF) of scalar fluctuations collapsed to 
the Gaussian distribution, whereas at large amplitudes it decayed more quickly than Gaussian. 
At large scales, the PDF of scalar increment behaved similarly to that of scalar fluctuation. In 
contrast, at small scales it resembled the PDF of scalar gradient. Further, the scalar dissipation 
occurring at large magnitudes was found to grow with Sc. Due to low molecular diffusivity, in the 
Sc > 1 flow the scalar field rolled up and got mixed sufficiently. However, in the Sc <g 1 flow 
the scalar field lost the small-scale structures by high molecular diffusivity, and retained only the 
large-scale, cloudlike structures. The spectral analysis found that the spectral densities of scalar 
advection and dissipation in both Sc 1 and Sc< 1 flows probably followed the fc -5 / 3 scaling. 
This indicated that in compressible turbulence the processes of advection and dissipation except 
that of scalar-dilatation coupling might defer to the Kolmogorov picture. It then showed that at 
high wavenumbers, the magnitudes of spectral coherency in both Sc 1 and Sc <S 1 flows decayed 
faster than the theoretical prediction of fc” 2 / 3 for incompressible flows. Finally, the comparison 
with incompressible results showed that the scalar in compressible turbulence with Sc = 1 lacked a 
conspicuous bump structure in its spectrum, but was more intermittent in the dissipative range. 


1. INTRODUCTION 

Turbulent mixing is of importance in many fields in¬ 
cluding the scattering of interstellar materials through¬ 
out the Universe, the dispersion of air pollutants in the 
atmosphere, and the combustion of chemical reactions 
within an engine [1-6]. In the literature of fluid dynam¬ 
ics, mixing in turbulent flows is often called as scalar 
turbulence. The related classical picture of cascade is 
that scalar fluctuations are generated at large scales and 
transported through successive breakdowns into smaller 
scales, the process proceeds until the scalar fluctuations 
are homogenized and dissipated by molecular diffusion 
at the smallest scale. Therefore, how a scalar gets mixed 
by a flow depends on whether its molecular diffusivity is 
small or large, even if the flow is fully turbulent. The 
common measure of diffusivity is based on the Schmidt 
number Sc = v/Xi where v and \ are the kinematic vis¬ 
cosity and molecular diffusivity, respectively. Generally, 
there are three different Sc regimes in turbulent mixing. 
For Sc ~ 1, previous experiments and simulations [7-9] 
have suggested that the scalar spectrum in the inertial- 
convective range follows the fc -5 / 3 power law when the 
Reynolds number is sufficiently high. The weakly diffu- 
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sive regime (Sc 1) has also received great deal of atten¬ 
tion [10-12], especially concerning the fc” 1 power law for 
the spectral roll-off in the viscous-convective range [13]. 
In terms of the strongly diffusive regime (Sc <C 1), re¬ 
cent simulations [14,15] have provided strong support for 
the putative theory proposed by Batchelor, Howells and 
Townsend (hereafter referred to as BHT) [16], namely, 
that the scalar spectrum in the inertial-diffusive range 
obeys the fc” 1 '/ 3 power law. 

As is well known, compressible turbulence is crucial to 
a large number of industrial applications and natural phe¬ 
nomena, such as the design of transonic and hypersonic 
aircrafts, inter-planet space exploration, solar winds, and 
star-forming clouds in a galaxy. Nevertheless, our current 
understanding of scalar transport in compressible turbu¬ 
lent flows lags far behind the knowledge accumulated on 
the incompressible one. Previous simulations of mixing 
in compressible turbulence using the piecewise-parabolic 
method [17,18] showed that for velocity, the compressive 
component is less efficient in enhancing mixing than the 
solenoidal component. Moreover, the scaling of scalar 
structure function accords well with the SL94 model. In 
this paper, we carried out a series of numerical simu¬ 
lations for compressible turbulent mixing, using a novel 
computational approach [19]. To examine in detail the 
effects of the Schmidt number on the scalar transport 
in compressible turbulence, the turbulent Mach number 
was fixed at around 0.30, whereas the Taylor microscale 
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Reynolds and Schmidt numbers were varied from 34 to 
216 and 1/25 to 25, respectively. For focusing the influ¬ 
ences brought by shock waves, a large-scale forcing with 
overwhelming compressive component is added to drive 
and maintain the velocity field [20]. This paper is part of 
a systemic investigation of the effects of basic parameters 
on compressible turbulent mixing. In a companion paper 
[ 21 ], we carefully examined the effects caused by changes 
in the Mach number and forcing scheme. 

The rest of this paper is organized as follows. The gov¬ 
erning equations and simulated parameter, along with 
the computational method used, are presented in Sec. 2. 
In Sec. 3, we first analyze the spectrum and structure 
function, then describe the probability distribution func¬ 
tion, and finally discuss the scalar transport in Sc 1 
and Sc <C 1 flows. The summary and conclusions regard¬ 
ing this paper are given in Sec. 4. 


2. GOVERNING EQUATIONS AND 
SIMULATION PARAMETERS 


We consider a statistically stationary system of a pas¬ 
sive scalar advected by the compressible turbulence of 
an idea gas. The velocity and scalar fields are driven and 
maintained by large-scale velocity and scalar forcings, re¬ 
spectively, where in the former the ratio of compressive 
to solenoidal components for each wavenumber is 20 : 1 . 
Furthermore, the accumulated internal energy at small 
scales is removed by cooling function at large scales. By 
introducing the basic scales of L for length, p 0 for den¬ 
sity, U for velocity, T 0 for temperature and <f>Q for scalar, 
we obtain the dimensionless form of governing equations, 
plus the dimensionless state equation of ideal gas, as fol¬ 
lows 
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V = pT. ( 2 . 5 ) 

The primary variables are the density p , velocity vector 
u, pressure p, temperature T and scalar </. The nondi- 
mensional parameters a and /3 are a = PrRe {7 — 1 )M 2 


and f3 = ScRe( 7 — 1 ) 7 . Tj is the dimensionless large-scale 
velocity forcing 
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(k /) exp(ik/x) + c.c., ( 2 . 6 ) 

1=1 

where Tj is the Fourier amplitude, it has a solenoidal 
component perpendicular to k; and a compressive com¬ 
ponent parallel to k;, and in magnitude the latter is 
twenty times over the former. Similarly, the dimension¬ 
less large-scale scalar forcing S is written as 

2 

S = Es(k*)exp(ik,x) + c.c.. ( 2 . 7 ) 

i -1 


However, there is only a solenoidal component in S. The 
details of the thermal cooling function A can be found 
in [19]. The viscous stress ay,- and total energy per unit 
volume E are defined by 

/ diii duj \ 2 / a\ 

(2 ' 8) 


E = 


(7-I ) 7 M 2 2 






(2-9) 


Here, 6 = duk/dxk is the velocity divergence or dilata¬ 
tion. M = U /Co is the reference Mach number and 
c 0 = y/jRT 0 is the reference sound speed, where R is the 
specific gas constant, and 7 = C p /C v is the ratio of spe¬ 
cific heat at constant pressure C p to that at constant vol¬ 
ume C v . We shall assume that both specific heats are in¬ 
dependent of temperature, which is a reasonable assump¬ 
tion for the air temperature in the simulation of current 
Mach number [19]. By adding the reference dynamical 
viscosity po, thermal conductivity kq and molecular dif- 
fusivity xoj we obtain three additional governing parame¬ 
ters: the reference Reynolds number Re = poUL/po , the 
reference Prandtl number Pr = poC p /no , and the refer¬ 
ence Schmidt number Sc = v 0 /xcb where vq = Po/Po is 
the reference kinematic viscosity. In current study the 
values of 7 and Pr are set as 1.4 and 0.7, respectively. 
Thus, there remain three independent parameters of M, 
Re and Sc to govern the system. For completion, we 
employ the Sutherland law to specify the temperature- 
dependent dynamical viscosity, thermal conductivity and 
molecular diffusivity as follows 


_ 1.4042T 1 - 5 
^ ,K,X ~ T + 0.4042' 


( 2 . 10 ) 


The system is solved numerically in a cubic box with 
periodic boundary conditions, by adopting a new com¬ 
putational method. This method utilizes a seventh-order 
weighted essentially non-oscillatory (WENO) scheme [22] 
for shock regions and an eighth-order compact central fi¬ 
nite difference (CCFD) scheme [23] for smooth regions 
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TABLE I: Flow statistics in the simulations. 


Case 

Sc 

M t 

Rex 

V 

Vb 

Voc 

v! 

</' 

Ek 

Ecf, 

<£> 

(e^) 

7> 

Cl 

25 

0.28 

35 

0.027 

0.005 

X 

2.06 

2.75 

2.13 

3.81 

1.20 

0.80 

0.37 

C2 

5 

0.28 

35 

0.027 

0.012 

X 

2.06 

2.44 

2.15 

3.00 

1.19 

0.87 

0.52 

C3 

1 

0.28 

34 

0.027 

0.027 

X 

2.06 

2.14 

2.11 

2.30 

1.21 

0.92 

0.70 

C4 

1 

0.29 

208 

0.007 

X 

0.007 

2.11 

2.24 

2.28 

2.53 

0.41 

0.90 

1.98 

C5 

1/5 

0.29 

207 

0.007 

X 

0.023 

2.11 

2.11 

2.24 

2.25 

0.40 

0.93 

2.31 

C6 

1/25 

0.30 

216 

0.007 

X 

0.078 

2.15 

1.98 

2.30 

1.97 

0.43 

1.25 

3.40 


outside shocks. A flux-based conservative formulation is 
implemented to optimize the treatment of interface be¬ 
tween the two regions and then improve the computa¬ 
tional efficiency. The details have been described in [19]. 
Instead of M and Re, the compressible flow is directly 
governed by the turbulent Mach number (M t ) and Tay¬ 
lor microscale Reynolds number (Re\) [24], which are 
defined as follows 


M t = M 


(Vry 


(2.11) 


Re\ = Re 


u'Mp) 

\/3(p) ’ 


(2.12) 


where v! = {uf) and A = v! / \J((duj / dxj) 2 ) are the 

root-mean-square (r.m.s.) velocity magnitude and the 
Taylor microscale, respectively. The sign (■) denotes en¬ 
semble average and the repetition on subscript stands for 
Einstein summation. Here, we point out that the defini¬ 
tions of M t and Re\ are based on nondimensional vari¬ 
ables, and this way is continuously used in the following 
text if there is no special illustration. 

Tabic [T] presents the major simulation parameters. The 
simulations are conducted on a N 3 = 512 3 grid and di¬ 
vided into two groups according to the value of Re\. The 
first group including Cl, C2 and C3 is used to study 
scalars with low molecular diffusivity in low Re flows, 
where Sc is decreased from 25, 5 to 1, and Re\ is around 
35. For the three cases, the smallest scale for velocity is 
the Kolmogorov scale p = [(p/(Rep)) 3 / < e/ p >] 1 ^ 4 , and 
that for scalar is the Batchelor scale ps = Sc~ l ^ 2 p with 
values of 0.005, 0.012 and 0.027. In the second group ad¬ 
dressing C4, C5 and C6, we pay our attention to scalars 
with high molecular diffusivity in high Re flows, where 
Sc is decreased from 1,1/5 and 1/25, and Re\ is around 
210. Here, the smallest scale for a scalar is the Corrsin 
scale, 77 c = Sc~ 3 ^ 4 p, instead of the Batchelor scale, pb, 
with values of 0.007, 0.023 and 0.078. 

The r.m.s. magnitude of velocity v! and the kinetic 
energy per unit volume Ek are increased by Re\ rather 
than by Sc. In contrast, the r.m.s. magnitude of scalar 
<p' and the scalar variance per unit volume E^ mainly 
grow with Sc in both high and low Re flows, where the 
related definitions are </>' = sj (c/ 2 ), Ek = (pu 2 ) / 2 and 



FIG. 1: (color online). Normalized scalar dissipation rate 
versus the product of Re\ and Sc 


E^ = (pc/ 2 )/2. Because M t is fixed, the ensemble-average 
value of the kinetic energy dissipation rate (e) becomes 
only dependent of Re In contrast, that of the scalar dis¬ 
sipation rate (e^) in both high and low Re flows increases 
as Sc decreases. Here, e = aijSij/Re, e^ = x(d(j)/dxj) 2 , 
and Sa = (dui/dxj +duj/dxi)/2 is the stain rate ten¬ 
sor. For the ratio of the mechanical to scalar timescales 

= (</OU 2 /2>/(e))/((p</> 2 /2)/<e0>), its dependence on Sc 
is similar to that of (e^). 

Further, an alternative to r $ is the quantity of //, = 
(Lf/u')/((p(j ) 2 / 2)/(cu)), where Lf is the integral length 
scale [25]. In Fig. fl] we plot against the logarithm 
of the product of Re\ and Sc. It shows that falls 
monotonously when the product increases. 


3. SIMULATION RESULTS 


A. Spectrum and Structure Function 


By applying the Kolmogorov theory [26,27] to the 
scalar transport in incompressible turbulent flows, 
Obukhov [28] and Corrsin [29] derived a scalar spectrum 
in the inertial-convective range satisfying L/) 1 « t < 


(3-i) 


E^(k) = C^(e 4 ,)(e)- 1 / 3 k- 5 / 3 , 
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FIG. 2: (color online). Compensated spectrum of scalar ac¬ 
cording to the Obukhov-Corrsin variables at different values 
of Sc and Re\. 



FIG. 3: (color online). One-dimensional compensated spec¬ 
trum of scalar at different values of Sc and Re\. 


where is the integral length scale of scalar [25]. C$ is 
the Obukhov-Corrsin (OC) constant, and the typical val¬ 
ues are 0.75 ~ 0.92 by experiments and 0.87±0.10 by sim¬ 
ulations [8,30]. In Fig. [2] we plot the compensated spec¬ 
tra of scalar according to the OC variables at different Sc 
and Re x . For the curves of E^k), plateaus appear in the 
inertial-convective ranges, especially for the Sc = 1 flows. 
This means that in the range of L7 1 « t « 77 1 , the 
scalar spectrum in compressible turbulent mixing seems 
to also obey the k ~ 5 / 3 power law. In high Sc flows, E c / ) (k) 
in the regime between the inertial-convective and dissipa¬ 
tive ranges exhibit as growing with wavenumber, which 
is reinforced by the increase in Sc. Contrarily, in low Sc 
flows, in the same regime E^k) falls as wavenumber in¬ 
creases, and this behavior is enhanced by the decrease in 
Sc. It implies that in a certain range, the scalar spectrum 
in a low or high Sc flow may have additional scaling. 

Although it is straightforward to compute in simula¬ 
tions the three-dimensional (3D) scalar spectrum as a 
function of wavenumber, experiments usually measure 
only the one-dimensional (ID) version of Ei^^k). In 



FIG. 4: (color online), (a) Compensated spectrum of scalar 
according to the Batchelor variables at Sc = 25, 5 and 1 and 
low Re where the slope value of the short line is —2/3. (b) 
The same as (a) according to the Batchelor-Howells-Townsend 
variables at Sc = 1, 1/5 and 1/25 and high Rex, where the 
slope value of line is —4. 


isotropic turbulence, it is written as 

OO 

Ei^k) = - J dk. ( 3 . 2 ) 

k 

Fig. [3] presents the ID compensated scalar spectra at dif¬ 
ferent Sc and Re\. We have taken averages over three 
coordinate directions. Previous studies of incompress¬ 
ible turbulence have showed the existence of a spectral 
bump, which is a precursor to the k~ l part of scalar 
spectrum, and becomes more and more pronounced as 
Sc increases. In our simulations, although it gets clearer 
when Sc grows, the bump is not as conspicuous as that 
observed in [31]. Similar to C^, the ID OC constant, 
Ci 0 , is changed by both Sc and Re\. We find that C 0 
and Citj, approximately satisfy the relation = 5 Ci 0 / 3 , 
which can be obtained directly through Eq. (3.2). 

When the Schmidt number is Sc 1, the energy spec¬ 
trum decays quickly at wavenumbers larger than r? -1 , 
whereas the scalar spectrum remains excited at levels 
higher than the energy spectrum [32]. In this case, the 
scalar transfer to small scales is creased at the Batch¬ 
elor scale, t]b , through molecular diffusion. The range 
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ry -1 CfcC rjg 1 is called as the viscous-convective range, 
wherein the scalar spectrum obeys a k^ 1 power law as 
follows 


= J B 0 (e 0 )(i//(e)) 1/2 fc \ ( 3 . 3 ) 

Here, the nondimensional coefficient B^ is presumed to 
be universal [31,32], and the value is 3 ~ 6 [33,34]. The 
nonlocality of the scalar transfer in wavenumber space 
is essential for the generation of this viscous-convective 
range. Therefore, a sufficiently high Sc is required to 
observe the fc -1 power law. In Fig. 0 a ) we plot the 
compensated spectra of scalar according to the Batche¬ 
lor variables in low Re flows, as functions of krj. For the 
Sc = 25 flow, a plateau representing the fc -1 power law 
is observed in the range of 0.1 < kr] < 0.3 (gray region), 
and the crossover region occurs in 0.06 < kr] < 0.1. How¬ 
ever, the fc -1 power law disappears when Sc is reduced to 
5. This reveals that in the viscous-convective range, the 
scalar spectrum in the Sc 1 compressible turbulent 
mixing defers to the k scaling given by the Batchelor 
theory, which was previously developed for incompress¬ 
ible turbulence. 

On the other hand, when the Schmidt number is 
Sc < 1, the scalar fluctuations at scales smaller than 
the OC scale, r/oc = Sc -3 / 4 ??, decay strongly, and thus, 
the scalar spectrum rolls off steeper than the k ~ 5 / 3 power 
law. If the Reynolds number is sufficiently large, stan¬ 
dard arguments suggest that in the range of iiqq <C k -C 
? 7 -1 , a so-called inertial-diffusive range exists. The afore¬ 
mentioned BHT theory predicts that in this range the 
scalar spectrum has the form 

E*{k) = (C K / 3 )(e^(e) 2 / 3 x~ 3 k- 17 / 3 . (3.4) 

Fig.gb) shows the compensated spectra of scalar accord¬ 
ing to the BHT variables in high Re flows, as functions 
of krj. We observe that in the Sc = 1/25 flow, a narrow 
plateau representing the k~ 17 / 3 power law arises, where 
the scale range is roughly 0.3 < kr] < 0.4 (gray region), 
which corresponds to 3.4 < kr]oc < 4.5. The above re¬ 
sult demonstrates that in the inertial-diffusive range, the 
scalar spectrum in the Sc -C 1 compressible turbulent 
mixing follows the fc -17 / 3 scaling provided by the BHT 
theory. 

The second-order structure function of scalar incre¬ 
ment is defined by 


5 ^ 2 (r) = ((< 5 r </>) 2 ), (3-5) 

where S r (f> = </>(x + r) — </>(x) is the scalar increment. In 
Fig. [5] we plot S^r) normalized by the OC variables, 
as suggested in Eq. (3.1), as functions of the normal¬ 
ized separation distance r/rj. It shows that finite-width 
plateaus at large r/rj. The scaling constants computed 
by the plateaus in low Re flows are higher than those in 
high Re flows. When Sc decreases, the values of scaling 
constants in both high and low Re flows become smaller. 



FIG. 5: (color online). Obukhov-Corrsin scaling of second- 
order structure function of scalar at different values of Sc and 

Re x - 



FIG. 6 : (color online). Yaglom scaling of mixed third-order 
velocity-scalar structure function at different values of Sc and 
Rex. 


At sufficiently large scales, an asymptotic formulation ob¬ 
tained from incompressible turbulent mixing [31] gives 

((<5r</>) 2 ) _ 3 RexSc 1 / 2 r 2/3 

(e<y,)(e) _ 1 / 3 r 2 / 3 \/l5 ? 7 oc 


3 Re\ ,r\—2/3 
\/l5 V 


(3-6) 


Given that the scale range for a plateau is different in 
every case, in principle it explains the relative relations 
of the scaling constants in our simulations. 

The mixed third-order structure function, defined as 
Sm 3 (r ) = ( S r u(S r 4> ) 2 ), where 5 r u = u(x + r) — u(x) is 
the longitudinal velocity increment, plays a more funda¬ 
mental role in the similarity scaling. In incompressible 
turbulence, an exact result for r/B <C r <C was given 
by Yaglom [35] 

(5 r u(5 r (j)) 2 ) = — ^ < £0 > r. ( 3 . 7 ) 

Fig.[ 6 ]presents the minus of S m 3 (r) normalized by the Ya¬ 
glom variables. It is found that for each simulated flow, 
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FIG. 7: (color online). The one-point PDF of the normalized 
scalar fluctuations at different values of Sc and Rex, where 
the circles are for the Gaussian PDF. 

there appears a flat region with finite width. Further¬ 
more, in the limit of large scale, — S m 3 (r ) drops quickly 
and approaches zero, while at small scales, it behaves ap¬ 
proximately as r 2 according to the Taylor expansion. We 
now employ quantity C m to represent the compensated 
mixed third-order structure function 

C„ = -4=V (3-8) 

Our results show that in flat regions, the scaling con¬ 
stant C m is 1.64, 1.43, 1.15, 1.16, 1.04, and 0.97 from Cl 
through C6. Two aspects can thus be concluded: (1) the 
contribution from the variation in the Reynolds number 
to C m is negligible, while the compressible effect makes 
C m smaller than the 4/3 value from incompressible tur¬ 
bulence; and (2) C m has a tendency to increase with Sc, 
which is in agreement with the results shown in [31]. In 
the inset of Fig. [6] we plot C m as a function of Sc. 

B. Probability Distribution Function 

Fig- □ shows the one-point PDF of the normalized 
scalar fluctuations. At small amplitudes the PDFs col¬ 
lapse to the Gaussian distribution, whereas at large am¬ 
plitudes they decay more quickly than Gaussian and thus 
are known as sub-Gaussian. This feature corresponds to 
the passive scalar transport in 3D compressible and ID 
incompressible turbulent flows [32,36], which is due to 
the fact that in our simulations, the ratio of the integral 
length scale of scalar L^ to the computational domain Lq 
is 0.19 ~ 0.21, which prevents large scalar fluctuations. 

In Fig. [8] we plot the one-point PDF of the normal¬ 
ized scalar gradient, where (/ = \J((dfi/dxj) 2 ) is the 
r.m.s. magnitude of scalar gradient. Obviously, the con¬ 
vex PDF tails are much longer than those of Gaussian, 
indicating strong intermittency. In both high and low 
Re flows, the PDF tails on each side become broader as 
Sc increases. Moreover, in the Sc = 1 flows, the notable 
increase in Re\ leads the PDF tails to be significantly 



FIG. 8: (color online). The one-point PDF of the normalized 
scalar gradient, where the circles are for the Gaussian PDF. 




FIG. 9: (color online). The two-point PDF of the normalized 
scalar increment. ( Sc,r/ri ) lines: (25,1), solid; (5,1), dashed; 
(1,1), dotted; (25,256), long dashed; (5,256), dash-dotted; and 
(1,256), dash-dot-dotted. The circles are for the Gaussian 
PDF, and the arrows indicate the decreasing Sc. (a) low 
Reynolds number, (b) high Reynolds number. 


wide. These observations imply that the growth in the 
Reynolds and Schmidt numbers will reinforce the events 
of extreme scalar oscillations at small scales. 

As a further study, in Figs. [9][a) and[9](b) we plot the 
two-point PDFs of scalar increment against the normal¬ 
ized separation distances of r/rj = 1 and 256, where cr scf, 
is the standard deviation of 5(j). We find that in both 
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FIG. 10: (color online). The PDF of the normalized scalar 
dissipation rate at different values of Sc and Re 


high and low Re\ flows, the behaviors of the PDF tails 
at large and small scales are respectively similar to the 
one-point PDF tails shown in Figs. [7] and [ 8 j 

Many studies in the literature of incompressible turbu¬ 
lent flows suggest that, at small scales, the intermittency 
of scalar field is closely associated with its dissipation. 
Fig. 10 shows the PDFs of the normalized scalar dissipa¬ 
tion rate, where e'^ = y/ ((e^ — (e^)) 2 ) is the r.m.s. mag¬ 
nitude of £ 0 . Similar to that shown in Fig. [ 8 j it is found 
that an increase in the Schmidt number strengthens the 
scalar dissipation occurring at large magnitudes. 



0 1 2 3 4 5 6 

X 


FIG. 11: (color online). Two-dimensional contours of scalar 
field in (a) Cl and (b) C6, at 2 = n/2. 


C. Scalar Transport Analysis: high Sc versus low Sc 


A central problem in turbulent mixing is the transport 
of scalar fluctuations in the inertial-convective range. For 
5c>lorS , c<l, this process should be also related to 
the viscous-convective or inertial-diffusive range. Fig. |11| 
shows the two-dimensional contours of scalar fields in the 
z = ir/2 plane for the Sc = 25 and 1/25 flows. In the 
highest Sc flow the low molecular diffusivity leads the 
scalar field to roll up and sufficiently mix. Nevertheless, 
in the lowest Sc flow the scalar field loses the small-scale 
structures by the high molecular diffusivity and leaves 
the large-scale, cloudlike structures. Given that there 
are density fluctuations in compressible turbulence, we 
introduce the density weighted scalar <F = y/pcj). The 
governing equation of scalar variance is then obtained by 
Eqs.(2.1) and (2.4) as follows 


d 

dt 






/3 ' dxj ’ 


Therefore, in a statistically stationary state, we only fo¬ 
cus the first three terms. 

In Fig. [12] we present the two-dimensional contours of 
the advection terms in the z = it/2 plane for the Sc = 25 
and 1/25 flows. In the highest Sc flow, the small-scale 
regions of extreme scalar advection distribute approxi¬ 
mately randomly in space, and display as thin stream¬ 
ers. For the lowest Sc flow, the small-scale structures 
are basically smeared by the high molecular diffusivity, 
and thus, the remaining structures are large in scale and 
exhibited as clouds. Fig. 13 shows the contours of the 
scalar-dilatation terms in the 2 = 7t/2 plane for the same 
flows. Because of the strong degree of compressibility in¬ 
duced by forcing, there appear large-scale shock waves in 
the scalar-dilatation contours. Besides, in the vicinity of 
a shock front, the scalar undergoes drastic changes. To 
observe the detailed structures at both small and large 
amplitudes, in Fig. 14 we plot the logarithms of the dis¬ 
sipation term in the z = t: / 2 plane, where the color scale 
is determined as follows 


X 1 d 2 

+ p y/pdx 2 ^ 2^' (3 ' 9 ^ 

Here, the terms on the right-hand-side of Eq. (3.9) suc¬ 
cessively stand for the advection, scalar-dilatation, dissi¬ 
pation and diffusion. In statistically homogeneous turbu¬ 
lence, the global average on the diffusion makes it vanish. 


V’ = logi 0 {D/D 1 ). (3.10) 

Here, D = x(< 9 </>/ dxj) 2 / '/3, and D' = ^ (( D — ( D )) 2 ) 
is the r.m.s. magnitude of D. The color changes from 
blue to red when dissipation increases. In the highest 
Sc flow, the contour shows that the extreme dissipation 
regions are sufficiently mixed and randomly distributed. 






















FIG. 12: (color online). Two-dimensional contours of the 
advection term of the scalar variance equation in (a) Cl and 
(b) C6, at z = 7t/2. 


In contrast, the high molecular diffusivity in the lowest 
Sc flow leads the contour to retain only the large-scale 
cloudlike structures. 

We now take attention to the spectral analysis of the 
transport of scalar fluctuations. First, Eq.(3.9) in Fourier 
space is written as follows 

$*(k) JU(k) = -$*(k)uW 7 $(k) - $*(k) 0 $(k) 

+ §<T(k)^V 2 <I>(k), (3.11) 

p Vp 

where the carets denote the Fourier coefficients, and the 
asterisks denote complex conjugates. In Fig. [15] we de¬ 
pict the log-log plots of the spectral densities of the ad¬ 
vection and dissipation terms from Eq.(3.11). It is ob¬ 
served that in Sc = 25 and 1/25 flows, there appear k~ 5 / 3 
power laws for both advection and dissipation. This in¬ 
dicates that although in compressible turbulent mixing, 
the presence of large-scale shock waves significantly af¬ 
fects the transport of scalar fluctuations, the processes 
of advection and dissipation may follow the Kolmogorov 
picture. Note that the process of scalar-dilatation cou¬ 
pling may not obey the Kolmogorov picture. The in¬ 
sets show the spectra of advection and scalar-dilatation 


FIG. 13: (color online). Two-dimensional contours of the 
scalar-dilatation term of the scalar variance equation in (a) 
Cl and (b) C6, at z = 7 t/2 . 


normalized by the dissipation spectrum. In the highest 
Sc flow, both the advection and scalar-dilatation spectra 
fall quickly as wavenumber increases. In contrast, in the 
lowest Sc flow the advection spectrum increases at large 
wavenumbers; however, the scalar-dilatation spectrum is 
positive at small wavenumbers but becomes negative at 
large wavenumbers, indicating mutually opposing scalar 
transfer processes. 

To determine the scalar spectra in both highest and 
lowest Sc flows, we compute the compensated scalar 
spectra relative to the kinetic energy spectrum. Ac¬ 
cording to the aforementioned theories, the results are 
fc -2 at Sc = 25 and k^E^/E^ at Sc = 1/25. It 
shows that plateaus appear for k~ 2 ^ 3 E ( j > /E k in the range 
of 38 < k < 51 (gray region with dashed edge lines) and 
k 4 E^/E k in the range of 23 < k < 34 (gray region with 
dash-dotted edge lines), respectively, which indicates 
that the related spectra dehned in the viscous-convective 
and inertial-diffusive ranges have flat regions located at 
relative large and small wavenumbers, respectively. Un¬ 
doubtedly, the kinetic energy spectrum for both flows 
have the inertial range of E k [k) = C k (e) 2 ^ 3 k ~ 5 ^ 3 [21], 
where C k is the Kolmogorov constant. This in turn yields 
E^ as follows 


E<j, oc Cfc(e) 2 / 3 fc _1 , Sc = 25; 


(3-12) 
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FIG. 14: (color online). Two-dimensional contours of loga¬ 
rithm of the dissipation term of the scalar variance equation 
in (a) Cl and (b) C6, at z = 7 t/2 . 


<x C fe (e) 2 / 3 fc- 17 / 3 , 5c =1/25. ( 3 . 13 ) 

It is useful to work with the spectral coherency defining 
as 

p u <t,{k) = E u ^(k)/{Eu^E^k)} 112 . ( 3 . 14 ) 

Here, E ui p(k) is the cospectrum of velocity and scalar, 
which is defined by [37] 

E u $(k) = J dS k (u( k) 0 *(k)>, ( 3 . 15 ) 


where the integral f dS k is taken over a spherical shell in 
wavenumber space. Lumley [38] proposed a k~ 11,3 power 
law for E u( f,(k) in the inertial-convective range, which 
was in good agreement with a recent simulation study 
[37]. Since both E k (k) and E^(k) scale as k~ 5 / 3 under 
similar conditions, the corresponding result for p u< j,{k) 
must be k ~ 2 / 3 . In Fig. 17 at high wavenumbers, the 
spectral coherency in the Sc = 25 and 1/25 flows decay 
faster than fc -2 / 3 . This deviation is mainly caused by the 
contribution from the coupling of scalar and dissipation, 
which provides a different power law. Contrary to the 
observations from incompressible turbulent mixing [15], 
here when Sc increases, the spectral coherency falls more 
slowly with wavenumber. Furthermore, the scattering of 
the spectral coherency points in the higher Sc flow is 
because of the lower Reynolds number. 


FIG. 15: (color online). Spectral densities of terms in the 
scalar variance equation in Fourier space, where the solid and 
dotted lines are for the advection and dissipation terms, re¬ 
spectively. Inset: advection (solid line) and scalar-dilatation 
(dashed line) spectra normalized by the dissipation spectrum, 
(a) Cl, (b) C6. 



FIG. 16: (color online). The relative compensated scalar spec¬ 
tra of fc -2 / 3 E^/Ek and k 4 E<f,/Ek in Cl and C6, respectively. 


D. Comparison Between Compressible and 
Incompressible Results 

In the final subsection we discuss the comparisons of 
certain results between compressible and incompressible 
turbulence. Herein, we focus solely on cases of Schmidt 
number at unity. In Fig. [18] we plot the compensated 
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FIG. 17: (color online). Spectral coherency of scalars in Cl 
and C 6 . The slope value of line is —2/3. 



FIG. 18: (color online). Compensated spectrum of scalar 
according to the Obukhov-Corrsin variables. The data from 
C4, [31] and [32] are denoted by the solid, dashed and dotted 
lines, respectively. 

scalar spectra according to the OC variables from C4, 
[31] and [32], where the values of Re\ are 208, 240 and 
258, respectively. For C4 the plateau centered at around 
kr) ~ 0.11 gives that C$ ~ 1.12. In contrast, the plateaus 
from the two incompressible flows are centered at lower 
wavenumbers, and the values of are 0.67 ~ 0.68. The 
spectral bumps appearing in the two incompressible flows 
are quite conspicuous. However, in the compressible flow 
it is difficult to identify a clear bump. Furthermore, ow¬ 
ing to the contribution of dissipation from shock waves, 
in the dissipative range, the scalar spectrum of compress¬ 
ible flow decays more quickly than its two incompressible 
counterparts. 

Another discriminative issue involves the scalar inter- 
mittency in dissipative range. A common used method 
for quantifying this intermittency is to compute the so- 
called intermittency parameter, /x^, through the auto cor¬ 
relation of scalar dissipation rate, namely, 

(e^,(x)e^(x + r)) ~ r~ M *. ( 3 . 16 ) 

Fig. [19] presents a log-log plot of the auto correlations 
of from C4 and [30], as functions of the normalized 
separation distance r/r]. It shows that the values of /x^ 


FIG. 19: (color online). Auto correlation of scalar dissipation 
rate, as a function of r/r], where the squares and circles are 
for the data from C4 and [30], respectively. 

from the compressible and incompressible flows are 0.73 
and 0.61, respectively. This means that compared with 
its incompressible counterpart, the scalar dissipation field 
in the compressible turbulent mixing is more intermittent 
by virtue of the contribution from shock waves. 


4. SUMMARY AND CONCLUSIONS 

In this paper, we systematically studied the effects of 
the Schmidt number on passive scalar transport in com¬ 
pressible turbulence. The simulations were solved numer¬ 
ically by adopting a hybrid approach of a seventh-order 
WENO scheme for shock regions, and an eighth-order 
CCFD scheme for smooth regions outside shocks. Large- 
scale predominant compressive forcing was added to the 
velocity field for reaching and maintaining a statistically 
stationary state. The simulated flows were divided into 
two groups. One was used to explore the scalar with low 
molecular diffusivity in low Re flows, where the Schmidt 
number Sc was decreased from 25 to 1, and the Taylor 
mircoscale Reynolds number Re\ was around 35. The 
other was addressed the scalar with high molecular dif¬ 
fusivity in high Re flows, where Sc was decreased from 
1 to 1/25, and Re\ was around 210. Our results showed 
that in both groups the ratio of the mechanical to scalar 
timescales increases as Sc decreases. As an alternative to 
r 4 >i f<p was found to fall monotonously when the product 
of Re x and Sc increases. 

In the inertial-convective range of L]/ 1 <C k <C rj~ l , 

the scalar spectrum seems to obey the k~ 5 / 3 power law, 
especially for the Sc = 1 flows. Besides, the tendencies of 
the scalar spectra to grow and fall with wavenumbers be¬ 
tween the inertial-convective and dissipative ranges are 
reinforced by the increase and decrease in Sc, respec¬ 
tively. For the ID scalar spectrum, the spectral bump 
becomes more visible as Sc increases, and the related 
OC constant and its 3D counterpart satisfy the rela¬ 
tion C 10 = 3C/i/5. Further, the scalar spectrum in the 
viscous-convective range of i; -1 C i: C % from the 














11 


Sc = 25 flow follows the k 1 2 power law, while that in 
the inertial-diffusive range of 77 /]^, <C k <C ? 7 _1 from 
the Sc = 1/25 flow shows a fc -17 / 3 4 5 6 7 8 9 10 11 scaling. At large 
scales, the scaling constant computed by the second-order 
structure function of scalar increment can be approxi¬ 
mately described using an asymptotic formulation devel¬ 
oped for incompressible turbulence. Simultaneously, the 
one computed by the mixed third-order structure func¬ 
tion of velocity-scalar increment shows that the depen¬ 
dence on the Reynolds number is negligible, and the ef¬ 
fect of compressibility makes it smaller than the 4/3 value 
from incompressible turbulence. In addition, this scaling 
constant has an increasing tendency when Sc grows. 

At small amplitudes, the one-point PDF of scalar fluc¬ 
tuations collapses to the Gaussian distribution, whereas 
at large amplitudes it is sub-Gaussian, exhibiting as de¬ 
cay more quickly than Gaussian. For the one-point PDF 
of scalar gradient, the convex PDF tails are much longer 
than Gaussian, indicating strong intermittency. In both 
high and low Re flows, the PDF tails on each side become 
broader as Sc increases. Furthermore, for the Sc = 1 
flows a notable increase in Re\ leads the PDF tails to 
be significantly wide, implying that the growth in the 
Reynolds and Schmidt numbers will enhance the events 
of extreme scalar oscillations at small scales. At large 
scales, the behavior of the two-point PDF of scalar in¬ 
crement resembles the one-point PDF of scalar fluctua¬ 
tion, while at small scales, it is similar to that of scalar 
gradient. In terms of scalar dissipation, it occurs more 
at large magnitudes when Sc grows. 

The contour shows that in the highest Sc flow the 
scalar field rolls up and gets sufficiently mixed, whereas 
in the lowest Sc flow it loses the small-scale structures be¬ 
cause of high molecular diffusivity, and leaves the large- 
scale, cloudlike structures. A further study on the con¬ 
tours of scalar advection and dissipation finds that in the 
highest Sc flow, streamers with extreme values are small 
in scale and are distributed approximately randomly in 
space, whereas in the lowest Sc flow only the large-scale, 
cloudlike structures exist. In certain ranges, the spectral 


densities of scalar advection and dissipation seem to have 
the A ; -5 / 3 power law. This indicates that for the trans¬ 
port of scalar fluctuations in compressible turbulent mix¬ 
ing, the advection and dissipation other than the scalar- 
dilatation coupling may follow the Kolmogorov picture. 
By computing the compensated spectra of scalar relative 
to the kinetic energy spectrum, it is confirmed that the 
scalings of k~ x and Ac 17 / 3 are defined for the scalar spec¬ 
tra in the viscous-convective and inertial-diffusive ranges, 
respectively. It then shows that at high wavenumbers, 
the magnitudes of spectral coherency in both the highest 
and lowest Sc flows decay faster than fc -2 / 3 , which is not 
similar to the prediction from classical theory. Finally, 
the comparison with incompressible results displays that 
the scalar in the Sc = 1 compressible flow lacks a conspic¬ 
uous bump structure in its spectrum near the dissipative 
range; however, it is more intermittent in this range. 

In summary, the above findings reveal that the change 
in the Schmidt number has pronounced influence on the 
small-scale statistics and field structure of passive scalar 
in compressible turbulence. Besides, although the tur¬ 
bulent Mach number used in current study is not very 
high, the effect of compressibility still affects scalar mix¬ 
ing in certain respects. A deeper investigation on this 
topic under much higher Reynolds and Schmidt numbers 
will be carried out in the near future. Also, we will ex¬ 
amine the effects of Mach number and forcing scheme on 
compressible turbulent mixing. 
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